Finger patterns produced by thermomagnetic instability in superconductors 
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A linear analysis of thermal diffusion and Maxwell equations is applied to study the thermo- 
magnetic instability in a type-II superconducting slab. It is shown that the instability can lead to 
formation of spatially nonuniform distributions of magnetic field and temperature. The distribu- 
tions acquire a finger structure with fingers perpendicular to the screening current direction. We 
derive the criterion for the instability, and estimate its build-up time and characteristic finger width. 
The fingering instability emerges when the background electric field is larger than a threshold field, 
E > E c , and the applied magnetic field exceeds a value -fffi ng oc 1/yE. Numerical simulations sup- 
port the analytical results, and allow to follow the development of the fingering instability beyond 
the linear regime. The fingering instability may be responsible for the nucleation of dendritic flux 
patterns observed in superconducting films using magneto-optical imaging. 



I. INTRODUCTION 

The thermomagnetic instability or flux jumping is 
commonly observed at low temperatures in type-II super- 
conductors with strong pinning. 1-4 The instability arises 
because of two fundamental reasons: (i) motion of mag- 
netic flux releases energy, and hence increases the lo- 
cal temperature; (ii) the temperature rise decreases flux 
pinning, and hence facilitates the flux motion. This 
positive feedback can result in thermal runaways and 
global flux redistributions jeopardizing superconduct- 
ing devices. This mechanism was understood in early 
works, 5,6 and later on the thermomagnetic instability was 
studied thoroughly (see Refs. 1-4 for a review). In par- 
ticular, the threshold magnetic field for the instability 
was calculated and its experimentally found dependence 
on temperature, sample dimensions, and the applied field 
ramping rate were explained. 

The conventional theory of the thermomagnetic 
instability 1,2 predicts "uniform" flux jumps, where the 
flux front is essentially flat. In other words, the spatial 
extension of the instability region tends to be maximal 
since small-scale perturbations are stabilized by thermal 
diffusion. This picture is true for many experimental con- 
ditions, however, not for all. Numerous magneto-optical 
studies have revealed that the thermomagnetic instabil- 
ity in thin superconducting samples results in dendritic 
flux patterns. 7-18 In course of the dendritic instability the 
flux forms narrow "fingers" of typical width 20-50 /im 
and length up to the size of the sample. Such a behavior 
clearly contradicts to the conventional theoretical con- 
cepts and needs elucidation. 

Few attempts to describe a nonuniform development of 
the thermomagnetic instability have been made. Among 
them is a numerical solution of thermal diffusion and 
Maxwell equations that can result in a rather nonuniform 
temperature distribution for a bulk superconductor. 1 
Molecular dynamics simulations of flux quanta motion in 



superconducting film 13 can model dendritic flux and tem- 
perature patterns similar to those found experimentally. 
However, these numerical results still lack analytical sup- 
port. In particular, it is still unclear what kind of spatial 
structure can be formed during a flux jump, and under 
what conditions. A similar problem was analyzed in a re- 
cent work 20 where the propagating flux front was shown 
to acquire a non-uniform spatial structure if its speed is 
higher than some critical value, and the conductivity is a 
strong function of flux density. In the present study it is 
shown that these assumptions are not necessary require- 
ments for a superconductor to develop nonuniform flux 
jumps. 



In the present paper the spatial pattern of the insta- 
bility in a bulk superconductor is studied using the con- 
ventional approach 1 ' 2 ' 5 - linear analysis of a set of dif- 
ferential equations describing small perturbations in the 
electric field E and temperature T. In contrast to the 
previous investigations, we allow the perturbations to 
vary in any direction, i. e., both parallel and perpen- 
dicular to the direction of the background current j and 
field E. In this way we determine the stability criteria 
and also estimate the instability build-up time. As a 
main result we find that the most unstable perturbations 
are in the form of narrow fingers perpendicular to the 
background field E and occur if E is larger than some 
threshold value. This shape prevents current adjustment 
across the perturbed region and, hence, yields the fastest 
perturbation growth. Too narrow fingers are, however, 
suppressed by the thermal diffusion. Thus, the typical 
finger size, \J n{dj c / dT)^ 1 / E, where n is the thermal 
conductivity and j c is the critical current density, is de- 
termined by the competition between the Joule heat jE 
and thermal diffusion, kV 2 T. 
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FIG. 1: Superconductor geometry. 



II. BASIC EQUATIONS 

We shall study the instability in the simplest geome- 
try i. e., in a superconducting slab placed in a parallel 
magnetic field, see Fig. 1. The slab fills the semispace 
x > 0, and the external magnetic field H is parallel to 
the z-axis so that the screening current j flows along the 
y-axis. The current and magnetic field distributions in 
the sample are determined by the Maxwell equation, with 
a proper boundary condition 

curlB = /i j, B| x=0 =/x H. (1) 

Here the local magnetic field in the flux penetrated part 
of the slab is assumed larger than the first critical field, 
and hence, to a good approximation B(x, y) = hqH(x, y). 
To find the temperature and electric field in a super- 
conductor the corresponding thermal and the second 
Maxwell equations should be used: 

C(dT/dt) = kV 2 T+jE, (2) 
curlE = -dB/dt, (3) 

where C is the specific heat. 

These equations should be supplemented by a current- 
voltage curve j = j(E, B, T). In type-II superconductors 
the j(E) dependence is strongly nonlinear. As a result, 
a quasi-static critical state with j w j c (B,T) is formed. 
This will be the initial state from which the instability 
evolves. For simplicity we use the Bean model, i. e. we 
neglect any B dependence of the critical current density 
j c . The exact form of the current- voltage curve, 

i=j(T,E)(E/E). (4) 

is not crucially important. The only important point 
is that the E(j) curve is very steep, and therefore its 
logarithmic derivative is large: 

Here a is the differential conductivity, 

a{E) = dj/dE. (6) 



At low electric fields, the E{j) curve is often approx- 
imated by a power law, i. e., n is assumed indepen- 
dent of E, and E oc j n . Our approach is applicable also 
to the flux flow regime at high electric fields. In that 
regime <j{E) = cry is the flux-flow Ohmic conductivity 
and n(E) = jc/cr/E cx 1/E. 

The key dimcnsionless parameter of the model is the 
ratio of thermal and magnetic diffusion coefficients: 

T = /ioCTK/C . (7) 

The smaller r is, the slower heat diffuses from the per- 
turbation region into the surrounding areas. Hence, one 
can expect that for smaller r; (i) the superconductor 
is more unstable, and (ii) the formation of instability- 
induced nonuniform structures is more favorable. 



III. PERTURBATION ANALYSIS 
A. Linearization of the problem 

We seek solutions of the equations presented above in 
the form, 

T + 5T(x,y,i), E + 5E(x,y,t), (8) 

where T and E are the background values. The back- 
ground field E may be created by ramping the external 
magnetic field H, or by other sources as discussed in 
Sec. VI. In practice E is nonuniform, but for simplic- 
ity we disregard its coordinate dependence. For a weak 
non-uniformity that can be justified using the method of 
Ref. 21, based on Wentzel-Kramcrs-Brillouin approxima- 
tion. In this approximation the non-uniformity results 
only in replacement some of local quantities by the ones 
averaged over x. Hence, we get only insignificant nu- 
merical corrections. Numerical simulations in Sec. V 
show that this conclusion also holds in the realistic situ- 
ation when the non-uniformity of E is induced by the by 
ramping the external magnetic field H. Similarly we ig- 
nore any coordinate dependence of the background tem- 
perature. This can be done if it satisfies the inequality 
T{x, y) — T <C T c — T, where T c is the critical tempera- 
ture of the superconductor, and T is the sample-averaged 
temperature before the instability build-up. 

From the symmetry of the problem E x = 0, while for 
the perturbation <5E both SE X and SE y in general do not 
vanish. The linearization of the E(j) in (4) yields 

<,_(g w+ ,„)! +A («_„-). m 

Since the vector E is parallel to the y-axis, one has in the 
linear approximation that SE = SE y , and as a result one 
finds 

S3={^ST + aSE y ^ +jc 5 -^. (10) 
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We shall seek perturbations in the usual form: 



ST = T*6exp(\t/t + ik v r) + ik x £). 



(11) 



SE, 



x,y 



Ee x , y exp(Xt/t + ik y rj + ik x £) , (12) 



where £ = x/w, n = yjw, and 



1 



~jc dT" 



aCT* 



{IqCFW , w 



CT* 

ml 



Here 9 and e x>y are the Fourier amplitudes, Re A is the di- 
mensionless instability increment, to is the characteristic 
time of adiabatic heating, which coincides with the mag- 
netic diffusion time for the length w, and w is the char- 
acteristic scale for the adiabatic instability. 1 The wave 
numbers k y and k x characterize the scale of the pertur- 
bation along the y and x axes, respectively. Since the 
sample is assumed infinite in the y direction, the k y is 
arbitrary, while k x is determined by the width of the flux 
penetrated region and the corresponding boundary con- 
ditions. 

Let us define the Fourier amplitude of the dimension- 
less current perturbation 5j/j c as i. Using Eqs. (10)— (12) 
one finds the components of the vector i in the form 



e x , i y = -9 + £ y /n. 



(13) 



Using Eq. (2) one obtains the equation for the temper- 
ature perturbation 9 as 



\9 = -T(k 2 y + k 2 x )9 + (i y + £ y )/n . 
We find from Eq. (14) 

e= (i + lAOfg 

nX + nr(k 2 + k 2 ) + 1 • 



(14) 



(15) 



Then, using Eqs. (1) and (3), we can rewrite the 
Maxwell equation for the perturbation as 



k x [k x e } = An i . 



(16) 



Using the relations (13) we cast Eq. (16) into the equa- 
tion set for dimcnsionlcss components of the electric field 
perturbation 



ky k x 
k 2 + An 



-k 2 x £ y + k y k x £ x = Xn(-6 + £ y /r 



(17) 
(18) 



Note that these equations together with Eqs. (13) 
and (15) provide continuity of the current perturbation, 
i. e., divi5j = 0, as required. Substituting Eqs. (15) 
and (17) in Eq. (18) one finds the following dispersion 
equation providing nontrivial solutions for £ y : 



l-\-T{kl+kl) _ fc| 

iX + nr(k 2 + k 2 ) + 1 k 2 + nX 



(19) 



The corresponding quadratic equation for X(k X7 k y ) has 
the form 



where 



A 2 + PA + Q = , 



P = k 2 x + kUn - 1 + r(k 2 y + k 2 x ) . 



n \ n x y n y 



The system is unstable if Re X(k x , k y ) > 0. 



(20) 



(21) 



B. Qualitative Analysis 

The dispersion equation becomes more transparent 
when the heat conductivity can be neglected, i. e. r = 0. 
Then, 



A 2 + \{kl + k 2 /n - 1) + (fc 2 - k 2 )/n = . 



(22) 



First, we notice that at k x = the system is always 
unstable. This is not surprising since such solutions cor- 
respond to the case of a sample with fixed transport cur- 
rent, i v = 0, heated by the electric field E under adia- 
batic conditions. In this case SE and ST grow with the 
maximal possible rate, A = 1, and the characteristic time 
of the instability build-up is to. 



T = 0.01 



T = 1 




FIG. 2: The contour plots for the instability increment 
Re \(k x , k y ) obtained from Eq. (20) for n = 10. The brightest 
areas correspond to the fastest growth of instability. For low 
r perturbations with a finite k y have the maximal increment, 
while for large r (strong heat diffusion), uniform perturba- 
tions with k y = would grow fastest. 



For a finite sample the k x is not arbitrary because of 
the boundary conditions at the edges of the flux pene- 
trated region. Only some particular k x satisfy the bound- 
ary conditions, which makes the system more stable. For 
example, for perturbations uniform in the y-direction 
(k y = 0) the instability develops only if k x < 1. However, 
if we set k y — > oo, then the system becomes unstable for 
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any k Xl and we again arrive at the maximal growth rate, 
A = 1. This result can be understood physically, if we 
take into account that infinite k y correspond to a per- 
turbation in the form of a narrow finger directed along 
the y axis, i. e., perpendicular to the current flow. The 
current flow remains unperturbed by an infinitcsimally 
narrow finger i. e., the condition i y = lest favorable for 
the stability holds, like for the case k x = 0. In the case of 
wider fingers, the current adjusts itself to the tempera- 
ture fluctuation, which slows down the instability growth. 
So, if one neglects the thermal diffusion, the narrowest 
possible fingers are the most favorable (k y — * oo), and 
the superconducting state is utterly unstable. 

The thermal diffusion evidently suppresses the insta- 
bility growth. The suppression is most effective for large 
k y . As a result, we obtain some optimal value of k v , for 
which the instability increment A is maximal. The exis- 
tence of such an optimal k y is evident from the contour 
plot of Re A calculated for t = 0.01, see Fig. 2, left. The 
dashed line shows k y providing the maximal Re A for a 
given k x . However, if r is larger then the heat diffusion 
fully dominates the instability development. In that case 
the maximal A corresponds to k y = 0, see Fig. 2, right. 

IV. RESULTS 

In this section we solve the problem more accurately, 
and start by establishing the proper boundary conditions. 



A. Boundary conditions 

From the above analysis it is clear that a finger struc- 
ture may appear only for r < 1. Consequently, we focus 
only on this case. Since the thermal diffusion is then 
slower than the magnetic diffusion, we can impose only 
the electrodynamic boundary conditions. This is equiv- 
alent to neglecting the heat flux in the x direction, i. c., 
the term rfc^ in Eq. (14) can be omitted. 

The magnetic field at the slab surface is equal to the 
applied field, hence the perturbation at the surface is 
zero, 5h z = at x — 0. The magnetic field has only 
^-component, thus from Eqs.(3) and (17) one obtains 
SE' oc Sh z , and the first boundary condition is 

SE' y = 0, x = . (23) 

This condition also means that the current does not flow 
across the sample surface, Sj x oc SE X = at x = 0. 

Let us specify the boundary conditions at the flux 
front, x = I. In the flux-free region, x > I, the elec- 
tric field decays on the scale of the London penetration 
depth, which is much smaller than any spatial scale of 
the problem. Therefore, the continuity of the tangential 
component of the electric field requires, 

5E y = 0, x = I. (24) 



These boundary conditions together with Eqs. (13), (14) 
and (16) are satisfied when SE y oc cos(k x x/w) with 

k x = (tt/2) (w/l) . 

Now we can search for solutions of Eq. (20) with this k x , 
and as before, when Re A > the system is unstable. 

B. Instability criterion and increment 

Let first consider the spatially uniform case where 
there exists a well-known criterion for the thcrmomag- 
netic stability. 1-6 With k y = and using t< 1. i. c. for 
very slow thermal diffusion we find from Eq. (20) that 
the system is unstable if k x < 1. For the Bean model, 
where I = H/j c , this is expressed as 

H > ff a diab = (tt/2) VCT*/h , (25) 

which is the commonly used adiabatic criterion for flux 
jumps. 



1 — i 




k x = 1 A ' ■ 




FIG. 3: The instability increment KeX(k y ) found from 
Eq. (20) for n = 10 and different k x . Top: fast heat diffusion, 
the maximal increment corresponds to uniform perturbations 
(ky = 0). Bottom: slow heat diffusion, the maximal Re A is 
found at a finite k y . 

Let us next consider cases of non-zero k y , and ana- 
lyze the behavior of Re X(k y ). Shown in Fig. 3 (top) are 
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plots for r = 1. For small applied magnetic fields the 
system is stable, see the curve for k x = 1 .1 . As the field 
increases, the flux penetration depth grows, and hence 
k x goes down. For k x = 0.7 the system becomes un- 
stable, i. e., solutions with Re A > arise. Note that 
the instability appears first at k y = 0. For higher fields 
{k x = 0.2), the instability range extends to large k y too, 
but the maximal Re A always corresponds to k y = 0. 22 
Therefore, for relatively large r the instability develops 
in a uniform mode. 

However, for smaller r the Re \{k y ) behaves differently, 
see Fig. 3 (bottom). The maximal Re A can here occur for 
a non-zero k y . Moreover, it is possible that the system 
is stable with respect to uniform perturbations, while 
unstable for perturbations with finite k y , see the curve 
for k x = 1.5 . This means that a non- uniform structure 
along the y-direction will be formed. 

When the applied magnetic field gradually increases 
from zero, the instability first starts for some particular 
k x = k* when Re A = only for one single value of k y = 
k*. This is the case for k x — k* = 2.1 in Fig. 3 (bottom). 
To find these k* and fc* one needs two conditions. The 
first one is 



Q(k*,k* v )=0, 



(26) 



which is a quadratic equation with respect to (k*) , and 
the second one is that the discriminant of this equation 
is zero. Using Eq. (21) and the fact that n > 1, we find 



1/4 



1 



K — \ ~ ) ~~i= > &x — 



1 



(27) 



The instability occurs at k x < k*, and for large n this 
instability criterion can be written as 



It follows directly that the border between the uni- 
form and non-uniform regimes is given by the criterion 
Rc \(k*,k y = 0) = 0. Using Eqs. (20) and (27) one can 
rewrite the criterion as r = 1/n. Rewriting this in di- 
mensional form we conclude that for 



E > E c = ^, Kj c /C 
the instability will evolve non-uniformly. 

V. SIMULATIONS 



(31) 




E > (tt 2 /4) (KT*/j c l 2 ) . 



(28) 



One can see from Fig. 3 (bottom) that the value of k y 
where Re A has the maximum depends only weakly on 
k x . Therefore, a good estimate for the finger width d y in 
the y direction is w/k*. Thus 



1/2 



1 



Edjc/dTj (2n)V4 



(29) 



Once we go from the instability threshold towards lower 
k x , the increment Re A quickly becomes of the order of 
unity. Thus, the characteristic time of the instability 
development is of the order of the adiabatic time, to- 
The aspect ratio of the perturbed region is 



k;/K^{2nf/\ 



(30) 



Note that it is independent of the thermal parameters, 
C, k,T*, and determined only by the shape of the E(j) 
curve. 

As was seen from Fig. 3, the instability will develop 
uniformly for r = 1, and non-uniformly for r = 0.01. 



FIG. 4: Evolution of the temperature T and electric field 
E y distributions produced by simulations that illustrate the 
formation of a finger structure. The instability was triggered 
by a uniform electric field, E y = Eo, switched on at t = 0. 
The images (a-d) correspond to the times t/to = 1.6, 3.0, 3.2, 
and 3.3, respectively. 

In order to visually illustrate the formation of non- 
uniform structures, and to verify the validity of the above 
analytical results, numerical simulations based on the 
Maxwell and thermal diffusion equations Eqs. (1), (2) 
and (3) were carried out. In the simulations we went 
beyond the linear approximation and considered the full 
non-linear E(j) curve, which was chosen to be 



E 



J 



°t + (ic/i)"" 1 jco/Eo 



(32) 



where j C Q and Eq are constants. This is one of possible 
smooth interpolations between the flux creep regime at 
small currents with E oc j n , and the Ohmic flux flow 
regime E = j/<Jf at high j. Here the flux flow resistivity 
is much higher than the characteristic resistivity in the 
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a b c d 




FIG. 5: Evolution of the temperature and flux density distri- 
butions produced by simulations. The instability is triggered 
by applying a magnetic field. The images (a-d) correspond to 
the times t/to = 45,46.1,46.15, and 46.25, respectively. The 
formation of fingers and their propagation into the flux free 
area is clearly seen. 



flux creep regime, a J 1 = 10 4 i?o/jc0- The temperature 
dependence of the critical current density is assumed to 
be linear, j c = j c0 [l — (T — T )/T*}. The electrodynamic 
boundary conditions are dE y /dx(x = 0) = (constant 
external magnetic field), E y (x = I) = 0, E x (x = 0, 1) = 0, 
with I = 2w, and the periodic boundary conditions in the 
y direction. Since the thermal boundary conditions are 
not of crucial importance at r <C 1, we used the simplest 
ones, T(x = 0) = To (ideal heat removal at the surface) 
and dT j dx(x = I) = (the symmetry condition in the 
middle of the slab). 

Our analytical results predict that the instability will 
form a finger structure if a sufficiently large uniform back- 
ground electric field is present. Therefore, for initial con- 
ditions we assume a uniform electric field, -E^O < x < 
Ut = 0) = Eq. To introduce some non-uniformity into 
the system small random values 5Tr <C T* were added to 
the initial temperature for every discrete node. The ini- 
tial temperature is then given as T(x, t = 0) = Tq + <5Tr. 
The key parameters r = 0.001 and n ~ h = 30 are 
specified at E = Eq. Their dependences on the electri- 
cal field are given by Eqs. (5), (6) and (32). Since now 
Eq = E c /{nT) 3> E c , the condition (31) is fulfilled, and 
the instability is expected to develop in a non-uniform 
fashion. This is indeed confirmed by the calculated evo- 
lution of E y and T distributions presented in Fig. 4. The 
numerical solution was performed on a grid of 140x70 
nodes using a simple-step integration method. 



One can see from Fig. 4(a) that at small times the 
distributions of E and T arc essentially uniform along 
the y axis. Then, a finger structure is emerging (b) with 
protrusions perpendicular to the electric field direction, 
as predicted by our previous linear analysis. The sim- 
ulations also show how this finger structure is evolving 
beyond the linear regime. We can see that the electrical 
field in some fingers grows faster so that relative differ- 
ence between the fingers increases (c). Eventually, the 
most intense finger takes over and dominates the entire 
E distribution (d). We believe that the reason for such 
behavior is the increase of the differential resistivity as E 
grows, Eq. (32). The growth of E is significant: the av- 
erage value E y — 1.7 Eq for (a), and 21Eq for (d). Note 
that this growth cannot be traced from the presented 
images only because the gray scale was optimized for 
each individual image to provide the best contrast. More 
detailed simulations showed that the instability growth 
slows down only when the increasing E reaches the in- 
flection point on the E(j) curve before entering the flux 
flow regime. 

Next, we carry out simulations with different initial 
and boundary conditions. We start from zero electric 
field, E(t = 0) = 0, and assume that a linearly in- 
creasing magnetic field is applied to the slab so that 
—dEy/dx = dH/dt = 0.03wj c /to at x = 0. The other pa- 
rameters are the same except that now the slab halfwidth 
is 6w. The right edge of the distributions shown in Fig. 5 
corresponds to the middle of the slab. One can see from 
Fig. 5(a) that for small H the flux penetrates in the con- 
ventional way, and a Bean-like profile of flux density is 
gradually advancing into the slab. When H and corre- 
spondingly E increase further, an instability sets in and 
leads to the formation of fingers (b-d). The finger struc- 
ture is apparent in both the B and T distributions, es- 
pecially on the later stages when only few intense fingers 
remain. Remarkably, the fingers tend to propagate into 
the flux free region, strongly distorting the flux front (d) . 
One can speculate that these growing fingers eventually 
will develop into a complex dendritic flux pattern ob- 
served by magneto-optical imaging. 7-18 

The instability criteria and its growth rate found from 
the simulations arc in a good agreement with our ana- 
lytical results. Moreover, the simulations demonstrate 
that the finger instability arises even if some assump- 
tions made in the derivation are relaxed. In particular, 
one does not necessarily need a strictly uniform back- 
ground E and T distributions as assumed in the deriva- 
tion. Furthermore, the background E and T distributions 
can also be non-stationary, which is always the case in 
a real experiment. In fact, in the simulations relevant 
to Fig. 5, where the instability was triggered by increas- 
ing the applied magnetic field, the E and T distributions 
were non-uniform and non-stationary. The formation of 
finger structure also turned out to be rather insensitive to 
the boundary conditions. We have also carried out simu- 
lations assuming that j c in Eq. (32) depends on the local 
B according to the Kim model, 23 j c (B) oc (B + |B|) _1 . 
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Uniform 



Fingering 




FIG. 6: Instability phase diagram in the plane magnetic field 
- electric field. The horizontal line corresponds to the adia- 
batic criterion for uniform jumps, Eq. (25). For E > E c , the 
instability has a finger structure, and the criterion is given by 
Eq. (33). 



With Bq = SfiowjciO) we found similar distributions, 
thus proving that the finger instability can arise also in 
cases with a B dependent E(j). 

The simulations presented here have some similarities 
with those by Aranson et al. 19 The main differences are 
that Aranson et al. started from a fully-penetrated state, 
the instability was nucleated by a local heat pulse, and 
j c was generally nonuniform. As a result, the obtained 
patterns of T distribution look different from ours. Nev- 
ertheless, they also found that the instability results in a 
non-uniform T distribution only at small r. 



VI. DISCUSSION 

The results obtained in this work can be graphically 
summarized by the instability "phase diagram" shown in 
Fig. 6. For small electric fields, E < E c , the conventional 
uniform instability is favorable, and the adiabatic insta- 
bility criterion Eq. (25) is applicable. For E > E c the 
fingering instability develops, with the instability crite- 
rion given by Eq. (28). Using the Bean model, H = jj, 
we obtain the finger instability criterion as 



H > H fing = (tt/2) V*T*j c /E , E > E c 



(33) 



Figure also shows the border between the regions of uni- 
form jumps and fingering instability for H > -ff a diab 
that was calculated from Eq. (20) using two conditions, 
dReX/dk y \ ky=k( =0, and ReX(k x ,k ) = ReX(k x ,0). 

Strictly speaking our analysis applies to the case r < 1 , 
which is equivalent to E > E c /n. For smaller electric 
fields a similar stability analysis can be made, taking into 
account the heat flux along the x axis. As expected, we 
found that for E < E c /n the uniform development of in- 
stability is always preferable. The instability criterion is 



given by the well-known dynamic criterion, that is highly 
sensitive to the external cooling conditions. 1 ' 2 However, 
in all cases the flux jump field decreases monotonously 
with E, as indicated schematically in Fig. 6. 

The finger instability occurs only at rather large back- 
ground electric field. This field can be created by differ- 
ent sources. For example, if the applied magnetic field in- 
creases with a rate H, an electric field E ~ HI ~ HH/ j c 
is generated. Thus, when increasing H with a constant 
rate we move in the phase diagram in Fig. 6 along a 
straight line staring from the origin. For small H, one 
crosses the instability boundary at E < E c , resulting in 
a uniform flux jump. For large H, the stability is de- 
stroyed for smaller H, and results in the formation of a 
non-uniform spatial structure. The predicted downturn 
of the H{E) instability line at large E can be checked 
experimentally 

Numerical estimates were made using typical parame- 
ters for low-temperature superconductors at helium tem- 
peratures: j c = 10 10 A/m 2 , C = 10 3 J/Km 3 , T* = 10 K, 
K = 10~ 2 W/Km, and n = 30. We then find the follow- 
ing values for the characteristic fields, -ff a diab ~ 0.1 T, 
and E c m 0.1 V/m, a finger width of d y m 3 fixn for 
E ~ E c , and a build-up time of the instability, to, in 
the //s range. These estimates are not far from those 
reported in experimental papers, namely dendritic fin- 
gers of width 20-50 /im, 8,11,14,16 and the instability build- 
up time of ~ 0.1 /is. 8 ~ 10 The criterion for the finger- 
ing instability E > E c can also be written down as 
a < <t c = C/nfiQK. Using the numbers above we find 



3 x 10 9 n 



9 O-l^-l 



, which is a reasonable value for 



the flux-flow conductivity. Correspondingly, for a = a c 
one obtains r = 1/n ~ 1/30. 

Note that the electric field, E Ci needed for the fin- 
ger instability to occur is not very small. As an esti- 
mate, the magnetic field ramp rate ^lqH that induces 
the electric field E c is of the order of 10 2 T/s for 
1=1 mm. Rates of similar magnitude are conventionally 
used for pulsed magnetization of superconducting per- 



manent magnets 



21. 25 



In experiments reporting the fin- 



gering instability 7-18 the ramp rates were much smaller. 
One should keep in mind however that the actual electric 
field can be much larger than it follows from the simple 
estimate [IqHI. The reason is a strong non- uniformity of 
the flux penetration both in space and in time, see for 
review Ref. 26. Hence, one can expect rather large local 
electric fields that last longer than the inverse instability 
increment, > 1 /is. Other sources of large electric fields 
include random fluctuations of the superconductor pa- 
rameters due to, e. g., relaxation of mechanical stresses. 
A very large electric field can be also created on purpose, 
e. g., by a laser pulse, which nucleates highly nonuniform 
flux distributions. 8-10 

In any case, it is rather difficult to meet the fingering 
instability criterion, E > E c . This might be the reason 
why fingering is hardly observed in bulk samples. We 
are aware of only one experimental work 7 where an indi- 
cation of the discussed fingering instability in relatively 
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thick samples (with thickness up to 2 mm) was obtained. 
Another possible reason for why such observations are 
few is that flux jumps in bulk superconductors are often 
complete, or almost complete. This means that the tem- 
perature rises close to T c in the entire sample, leading to a 
uniform flux distribution which erases any trace of a pos- 
sible non-uniformity in the first stages of the flux jump. 
That contrasts the behavior in thin film samples, where 
the jumps are usually much smaller and far from being 
complete. 27 This makes non- uniform jumps easier observ- 
able in films. 8-18 Moreover, huge stresses usually exist 
between a superconducting film and substrate. Abrupt 
relaxation of these stresses can lead to fluctuations in E, 
especially for NbsSn or MgB2 where the superconducting 
properties depend strongly on the strain. 2,29 Note also 
that in films it is much more probable that any pertur- 
bation of electric field will influence the whole thickness, 
whereas it will affect only a small part of a bulk sam- 
ple. Although our equations for a slab cannot be directly 
applied to the case of a thin film, we expect that essen- 
tially the same physics describes the formation of finger 
structures in the films too. 28 Moreover, non- locality of 
the current-field relations in films can make formation of 
non-uniform structures there even more favorable, and 
possibly account for the branching flux patterns observed 
exp er iment ally. 

The presence of a background electric field E, and 
hence moving magnetic flux implies that the background 
state itself is not stationary. In a typical experiment, the 
applied magnetic field is increasing H = H(t), the flux 
front is moving into the sample, I = l(t), and hence the 
electric field is non-stationary within the flux-penctratcd 



region. Obviously, our analytical results are valid only if 
all these quantities change in time slower than the pertur- 
bations 5E, 8T grow, i. e., when E/E, H/H, l/l <C A. If 
the electric field is created by ramping the external mag- 
netic field, E(x,t) ~ fJ, H[l(t) - x], then E/E w H/H w 
l/l. Using that A ~ 1/ioj we can rewrite the above in- 
equality as H 3> i?adiab/ ' y/n. Since n> 1 this condition 
is satisfied in the major part of the phase diagram in 
Fig. 6. 

In conclusion, a linear analysis of heat diffusion and 
Maxwell equations shows that a thermomagnetic insta- 
bility may result in finger-like distributions of T , E and 
B. The fingering instability arises if the background elec- 
tric field is so high that the magnetic flux diffusion pro- 
ceeds much faster than the heat diffusion. Numerical 
simulations have shown that upon further development 
of the instability one finger starts growing much faster 
than the others, and propagates into the flux-free region. 
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